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Abstract 

In previous works, in order to include correction by the Coulomb wave function 
in Bose-Einstein correlations (BEC), the two-body Coulomb scattering wave func¬ 
tions have been utilized in the formulation of three-body BEC. However, the three- 
body Coulomb scattering wave function, which satisfies approximately the three-body 
Coulomb scattering Schrodinger equation, cannot be written by the product of the 
two-body scattering wave functions. Therefore, we reformulate the three-body BEC, 
and reanalyze the data. A set of reasonable parameters is obtained. 

1 Introduction 

Recently, in addition to the data on the two-body charged Bose-Einstein correlationsiBEC), 
data on the three-body charged BEC have been reported PPI3]. In some papers PEI, 
the Coulomb correction is done with fixed source radius, for example, 5 fm. On the other 
hand, the quasi-corrected data (raw data with acceptance correction) on the two-body 
(27t“) BEC 3} and the three-body (37t - ) BEC have been reported [2|. 

In Ref. [SIP, authors proposed a theoretical formula for the 37 t _ BEC by the use of the 
two-body Coulomb wave functions, and outputted information on BEC with fixed source 
radii (5 fm and 10 fm). On the other hand, we have analyzed the 2ir~ and 37 t~BEC, using 
the CERN-MINUIT program with the two-body Coulomb wave functions and the source 
radius as a free parameter [hi Emm 

The formula for 27r“BEC reduces to that of plane wave formulation in the limit of 
plane wave approximation. However, the formula for 37r _ BEC does not reduce to that 
of plane wave formulation mm Additional factor (3/2) appears in the phase of plane 
wave [HI El EDI- Therefore, we have re-interpreted the source radius estimated from the 
analysis of 37r~BEC. 
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In this paper, we would like to examine the relation between the two-body Coulomb 
wave function and the asymptotic solution of the three-body Coulomb wave function, which 
cannot be written by the product of two-body Coulomb wave functions. In addition, we 
would like to show that factor 3/2 disappears from the phase factors of plane wave in the 
formulation of the 37 t“BEC, if the correct asymptotic three-body Coulomb wave function 
is used. 

In the second section, an asymptotic solution for the three-body Coulomb scattering 
Schodinger equation is shown. The formula for 37 t _ BEC is derived from the analogy of the 
formula for the plane wave formulation in the third section. Analysis of 37r _ BEC is done 
in the fourth section. Final section is devoted to summary and discussions. 


2 Approximate solution for Schodinger equation of three- 
body Coulomb scattering 

In order to describe the two-body charged BEC ( for example, 27 t~ system ), we should 
solve the Shrodinger equation of Coulomb scattering. The solution, which is regular at the 
origin of the Coulomb potential, is given by 

= < ' k:J ' ri ’Ok. , [r;/). 

<l>kij( r ij) = r ( 1 +iVij)e’ rT,ij/2 F[-irii j ,l;i(\ki j \-\r ij \-k ij -ri j )}, ( 1 ) 

for particles i and j, where the coordinate and momentum of particle i, are denoted by 
xi and fcj, respectively, and ej in is the charge of particle i. In Eq. (p~t . the relative 
coordinate and momentum of particles i and j are denoted by r t j = Xi — Xj. and kij = 
(■rrijki — rriikj)/{m,i + rrij ), respectively, rjij = eiejiiij/\kij\ where Hij is reduced mass of 
rrii and rrij, F[a, 6; x\ is the confluent hypergeometric function, and T(x) is the Gamma 
function. 

In order to describe the three-body Coulomb scattering, the Jacobi coordinates M are 
introduced; 

Ci = X 2 ~Xi, 

C 2 = x 3 - (mixi + m 2 x 2 )/M 2 , 

C 3 = {mixi + m 2 x 2 + m 3 x 3 )/M , (2) 

M 2 = mi + m 2 , M = mi + m 2 + m 3 . 


The relative coordinates are written as, 

r 2i = x 2 X\ = Ci j 

^31 = x 3 - Xi = aCi +C2, 

r 3 2 = x 3 -x 2 = -( 3 Ci + C 2 , 

a = m 2 /M 2 ,/3 = rrii/M 2 . 


The Schrodinger equation of the three-body Coulomb scattering is given by, 
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Figure 1: Jacobi coordinates of three-body system. The starting point of </ 2 is the center 
of mass of particles 1 and 2. 


where 

Hi = mim 2 /M 2 , I-I 2 = M 2 m 3 /M, 

Pi = HidCi/dt = (mifc 2 — m 2 Aq)/M 2 , 

Po = H 2 dC 2 /dt = {M 2 k 3 - m 3 (fci + k 2 ))/M. 


Then, the approximate solution for the Schrodinger equation in flo, where r 32 , r 23 , r 3 1 >> 1, 
is given by C11II3, 

% = ed p ^+ p ^)(j) kl2 (n 2 )0 fc23 (r 23 )^fc3i (r 31 ). (4) 

The phase factor of the plane wave in Eq. m is rewritten as, 


+ m 2 m 2 + m 3 , m 3 + mu 

P 1 C 1 +-P 2 C 2 = - 77 - &12 • ^12 H- 77 - k 23 ' ^23 H- 77 -«31 • r 31, 


Af 


M 


M 


= -{kv2-r 12 + k23-r 2 3 + k 31 -r 31 ), 


where rri\ = m 2 = m 3 is used. 

Therefore, the solution Ty is written as US, 

^7 = X 2)^k2k 3 ( x 2,X 3 )i^' 3kl (x 3 , Xl ) 

= e* (2/3)fcij <f>k i:j (r i: j ). 


(5) 

( 6 ) 


The approximate solution for the three-body Coulomb scattering can be written by the 
product of but not the product of the wave function of two-body scattering, 

V’fcifc-(* 1 , x j)- I n the correct formula, factor 2/3 is multiplied to the phase of plane wave. 

In Ref.jZ], the Coulomb correction for n-body scattering is discussed, where the n- 
body Coulomb scattering wave function is given by the product of two-body Coulomb 
scattering wave functions. However, the n-body Coulomb scattering wave function in flo 
is approximately given by 
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for n > 3. 


3 formula for third order BEC 

The wave function of identical Bose particles should be symmetrized. In FigJ^l or FigEl 
V c denotes the interaction between two particles by the Coulomb potential, and cross (X) 
represents the exchange of particles. 

As is shown in FigH the two particle momentum density is given by, 



2 

^ ] = \f\.J P( X i) d3x i\' l Pk 1 k 2 (xi, x 2 ) +'0fc 1 fc 2 (*2, *l)| 2 

2=1 



2 , 

= I p(x i )^x i (G 1 + G 2 ), 

2=1 J 

Gl = ^ (l^fc lfc2 (^l, *2)| 2 + |^fc lfc2 (*2, «l)| 2 ) , 



G 2 = Re(^fc lfe2 (*i, x 2 )^* k2 (x 2 , *i)) , 


where 
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Figure 2: Two-body BEC diagram 


The exchange diagram for 3ir BEC is shown in FigEI According to the diagram, the 
three particle density for 37r~ BEC is written as, 


where 
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1 
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p(xi)d 3 Xi\J2 A(J)\“ 


3 = 1 


A(l) = Ax = ip% h 2 (xi, *2)V’£k 3 ( a! 2, x 3 )^Z' 3 kl (x 3 , *i), 
A{ 2) = A 23 = ^kik 2 ( x i, x 3 )^k2k 3 ( x 3, tf 2 )V’£fci(*2, *i), 
A ( 3 ) = A 12 = V’fc 1 , fc a (*2, *i)^tefc 3 (*i» BaWZkito, *2), 


( 8 ) 
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,4(4) = A 123 =i’k 1 k 2 (X 2 , X 3 )^ k3 (x 3 , *i)^g fcl (a!i, * 2 ), 

41(5) = ^132 =1pZk 2 ( x 3, asOV’SfcsC*!’ x 2 )^Zki (* 2 , * 3 ), 

,4(6) = A 13 = V’feife,^, x 2 )^k 2 k 3 {x 2 , *i)^£ fcl (*i, * 3 )- ( 9 ) 
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Figure 3: Three-body BEC diagram 


In the plane wave approximation, each amplitude A(i) approaches to the following 
form, 

A(l) = Al e i(2/3)(fcl2-r-l2+fc23-7’23+fe31-T3l) _ gi(fel•£Cl+fe2■ £E 2 + &3' *3 ) 

y4(9) = ylog P^ e i(2/3)(fci2-r-i3+fe23-r-32+fc3l-r-2l) _ e i(ki xi+k2-X3+k3 X2) 

2 ^ 3 ^ _ e *(2/3)(fci2 r2i+fc23-ri3+fc3i r32) _ e *(fci £C2+fc2-a:i+fc3 3 : 3 ) 

y4(4) = _4i23 e *( 2 / 3 )( fc 12-7’23+fc23 'f’31+fc3l-ri2) _ e i(k\-X2+k2X3+k3-X\) 


y4(5) = ^ 132 P ^y e i(2/3)(fci2-r 3 i+fc23 ri2+fc3i-r23) _ e i(fci-X3+fc2 a:i+fc3-cc2) 

y4(6) = y4 13 P ^y e i(2/3)(fcl 2 r32+fc23-r21+fc3i ri 3 ) _ gi(fcl ■IT3+fc2-®2+fc3 , ®l) 


( 10 ) 


In Eq. EH, PW means the plane wave approximation of the amplitude, and the condition 
in the center of mass system, exp[— i(k\ + k 2 + £ 3 ) • £ 3 ] = 1 is used. 

The amplitudes squared in Eq.® can be classified into the following groups, 


Fl = — [Al A* + 2 I 122 I 32 + ^23^23 + ^.132133 + ,4123,4^23 + VI 1322 I 132 ], 

F\2 = —[^1^12 + ^23^123 + ^134li32 + C.C.], 

F'n = “[ 24 i 24 23 + Ai 2 A ^ 32 + ^13^1*23 + C.C.], 
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^31 = “ [-A 1 -A 13 + ^423^132 + ^12^123 + C.C.\, 

F\ 23 = g[^4l^4l32 + ^132^123 + ^13^12 + ^12^23 + ^23^13 + ^4l23^l]) 

Tj32 = ^[^1^123 + ^23^12 + ^12^13 + ^123^132 + ^132^1 + ^13^23]! (H) 

where, c.c. denotes the complex conjugate, and T 132 is the complex conjugate of T 123 . In 
the plane wave approximation, T\ reduces to 1, T/,- corresponds to exchange between i and 
j charged particles, and T 123 correspond to exchange among three charged particles. 

Phenomenologically, the coherence parameter A is introduced into the formula for 
27 t~BEC as, 

7 V 2it_ 2 r 

Wi m = C Uj P(xi)d 3 Xi(G 1 + AG 2 ), (12) 

i=l' 

where C is the normalization factor. 

In the third order BEC, factor A ”/ 2 is multiplied to the amplitudes squared according 
to the number n of exchange particles. After £3 integration, the 37r“BEC is given by 
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AT BG 
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F\ + 3 AEi 2 + A 3//2 Re (F 123 ) . 

(13) 


The set of following variables is used in the concrete calculations of Eq. jm 

k\2 = —Pi, 

&23 = \Pl - \ P 2, 

&31 = \P\ + \P%, 

Q, = ^{k{ 2 + k \ 3 + kl) = ^6P 2 + |P 2 . 


(14) 


4 Analysis of 3 tc BEC 

The formula m is applied to the analysis of quasi-corrected data on 37 r BEC by STAR 
Collaboration j^j. The results are shown in Tabled and Fig. d For comparison, the result 
of previous work [Sj is also shown in the lower part of Tabled The result for 27r~BEC by 
STAR Collaboration [3J is shown in Tabled 

The source radius R% r d estimated from the data on 37 r _ BEC with Eq. Cl is comparable 
with that i? 2 nd from 27 t“BEC. However, the source radius R^ d estimated in the previous 
work d, namely with the two-body Coulomb wave functions is much smaller than R^ r d 
with Eq. m- The re-interpreted radius becomes (3/2 )R^ c d = 8.01 [fm], which is nearly 
equal to Rs r d- The coherence parameter estimated from 37 t~BEC with Eq. m is 
somewhat smaller than that \ 2 r 1 d from 27 t~BEC. 
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Table 1: Analyses of 3ir BEC by STAR Collaboration (Jjl. 



R [fm] 

A 

X 2 /n.d.f. 

Ea.ljT3l) 

8.26T0.39 

0.50T0.02 

1.88/35 

previous work !) 

5.34T0.24 

0.56T0.02 

2.80/35 


Table 2: Analyses of 2ir BEC by STAR Collaboration |1'. 

R [f'rn] A x 2 /n.d.f. 

Eg.(fT2|) 8.75T0.31 0.58T0.02 23.0/25 



Figure 4: Analysis of quasi-corrected data on 3n BEC by STAR Collaboration [&] with 

Eq.tQJ. 


The problem on the phase factors appearing in the two-body BEC among three identical 
particles muni is proposed in HU. If these factors are taken into account, the formula 
for 37 t~BEC is given by, 


N 3ir 
-J^BG 



1=1 


p(xi)d 3 Xi 


F 1 + 3XF 12 + 2A 3/2 Re[Fi 23 ] x W 


(15) 


where W = cos(</>i 2 + </>23 + ^ 3 i), which is parameterized as W = cos (g x Q%) in the simplest 
form with parameter g and Q 2 = (k\ — k 2 ) 2 + (k 2 — k^) 2 + (fc 3 — k\) 2 . The result is shown 
in Table Inland Fig. 21 Estimated source radius shown in Table 01 is smaller than R% r d with 
Eq.tEJ, and is not consistent with R 2n d- 


Table 3: Analyses of 37 t BEC by STAR Collaboration with Ea. (11-51) . 

R [fm] A ~g y 2 /n.d.f. 

Eq.(fT5l) 7.70T0.57 0.55T0.05 31.13± 10.87 0.51/34~ 
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Figure 5: Analysis of quasi-corrected data on 37 t BEC by STAR Collaboration }5 with 
Ea. (TF)|) . 


5 Summary and discussions 

We reformulate the formula for 37r“BEC, using the asymptotic three-body Coulomb wave 
function. We apply the formula to the analysis of data on 37 t~BEC by STAR Collaboration. 
The source radius R?, r d estimated from 37 t~BEC is consistent with that from 27r”BEC 
m ■ The coherence parameter X^ r d estimated from 3ir BEC with Eq. m is almost the 
same with A 2 n d from 27r _ BEC. Whether a set of preferable parameters can be estimated 
from the analyses of 27 t _ BEC and 37 t~BEC or not in other approaches will be reported 
elsewhere m- 

By the use of our formula, we can estimate source radius with Coulomb correction, 
without the re-interpretation due to the factor (3/2). 
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